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ABSTRACT 

This paper describes the Chandra observation of the diffuse emission in the face-on spiral 
NGC 6946. Overlaid on optical and Ha images, the diffuse emission follows the spiral structure 
of the galaxy. An overlay on a 6 cm polarized radio intensity map confirms the phase offset of 
the polarized emission. We then extract and fit the spectrum of the unresolved emission with 
several spectral models. All model fits show a consistent continuum thermal temperature with a 
mean value of 0.25±0.03 keV. Additional degrees of freedom are required to obtain a good fit and 
any of several models satisfy that need; one model uses a second continuum component with a 
temperature of 0.70±0.10 keV. An abundance measure of Sj^J gQ for Si differs from the solar value 
at the 90% confidence level; the net diffuse spectrum shows the line lies above the instrumental 
Si feature. For Fe, the abundance measure of 0.67±0.13 is significant at 99%. Multiple gaussians 
also provide a good fit. Two of the fitted gaussians capture the O Vll and O Vlll emission; the 
fitted emission is consistent with an XMM-Newton RGS spectrum of diffuse gas in M81. The ratio 
of the two lines is <0.6-0.7 and suggests the possibility of non-equilibrium ionization conditions 
exist in the ISM of NGC 6946. An extrapolation of the point source luminosity distribution 
shows the diffuse component is not the sum of unresolved point sources; their contribution is at 
most 25%. 

Subject headings: galaxies: individual: NGC 6946; galaxies: spiral; galaxies: ISM; X-rays: galaxies 



1. Introduction 

The total X-ray emission from a spiral galaxy 
consists of the sum of the emission from its X- 
ray-emitting constituents including X-ray binaries 
(low-mass and high-mass), supernova remnants, 
OB star clusters, other point sources, and the hot 
component of the interstellar medium. This last 
component is expected to be present based upon 
the theoretical work of McKee & Ostriker (1977) 
and Norman & Ikeuchi (1989) among many oth- 
ers. Prior to the launch of Chandra and XMM- 
Newton, the hot component had been detected 
in a handful of galaxies, largely as a consequence 
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of the broad spatial resolution of previous X-ray 
telescopes which distributed X-rays from point 
sources over a wide region of the galaxy, contami- 
nating any signal from the hot ISM. 

The study of the X-ray emission from normal 
galaxies commenced with the Einstein observa- 
tions of galaxies (for a review, see Fabbiano 1989). 
The spatial resolution of Einstein, ~1' (e.g.. Hart- 
line 1979), prevented separation of the X-ray emit- 
ting constituents for all but the nearest galaxies. 
ROSAT improved the resolution by about three 
times, providing views of a broader range galaxies 
outside the Local Group (e.g., Roberts & Warwick 
2001). The use of Chandra and XMM-Newton has 
led to more detections because of the large effec- 
tive area of XMM-Newton and the sharp point 
spread function of Chandra. 

NGC 6946 is a well-studied face-on spiral ini- 
tially detected in X-rays with the Einstein Obser- 
vatory (Fabbiano & Trinchieri 1987). Schlegel 
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(1994) obtained a 36-kscc pointing using the 
ROSAT PSPC, detected 9 point sources, and re- 
ported the existence of a diffuse component with a 
temperature of 0.55 keV (90% error range of 0.3- 
0.65 keV) which we will show just overlaps with 
the results to be presented. That study argued 
that the diffuse component was truly diffuse and 
not the summed emission of point sources based 
upon a plausibility argument for the expected 
numbers of sources of various types. In a separate 
study, Schlegel, Blair, & Fesen (2000) described 
a ROSAT BRl observation as well as an observa- 
tion using ASCA. In this case, the high internal 
background of the HRI prevented a detection of 
the diffuse component, while the broad PSF of 
the ASCA mirrors blended events from the known 
point sources with the diffuse emission. 

This paper describes the diffuse emission in 
NGC 6946 revealed by a 6G-ksec Chandra obser- 
vation of the galaxy. The point source population 
has been described by Holt et al. (2003). The pri- 
mary goal of this paper is a spectrum of the diffuse 
emission. Following this introduction, we describe 
our analysis steps, then present a brief compari- 
son of the diffuse emission with images from other 
bandpasses, followed by a longer analysis of the 
spectral fits and a discussion of the results. 

Chandra presents a unique opportunity to ob- 
serve the diffuse component as free as may be pos- 
sible, at least for the foreseeable future, from con- 
tamination by point sources. The Chandra mir- 
rors deliver ^80% of the X-rays at 1 keV into a 
circle ~0.7 arcsec in radius and 70% into a 0.7 
arcsec radius circle at 6 keV. At 5' off-axis, the 
corresponding numbers are 90% into a circle of a 
radius ~2".5 at 1 keV and ~4" radius at 6 keV 
(CXC Proposers' Guide 2002). 

The adopted distance for NGC 6946 is 5.9 Mpc 
(Karachentsev et al. 2000) with a corresponding 
scale of 1" ~29 pc. The column density in the di- 
rection of NGC 6946 is -3-5x10^1 cm'^ (Schlegel, 
Finkbeiner, & Davis 1998). 

2. Observation Summary 

Chandra observed NGC 6946 on 2001 Septem- 
ber 7. The galaxy was imaged on the back- 
illuminated (BI) chip S3. The data were exam- 
ined for flares to which the BI chip is susceptible, 
but no flares with a count rate greater than 0.04 



counts s^^ were located; screening any flare- like 
spikes removed 1296 sec of the 59116 sec observa- 
tion, leaving a net on-source, deadtime-corrected, 
exposure of 57819 sec. 

The data were filtered to eliminate events with 
energies below 0.3 keV and above 8 keV. The low- 
energy cut does not remove source photons be- 
cause those events are lost to the high column in 
the direction of the galaxy. The high cut removes 
extraneous cosmic ray-generated events. 

While the internal background of the ACTS de- 
tector is low (for S3, ~0.3 counts chip~^, CXC 
Proposers' Guide (2002)) and often may be ig- 
nored for counts from point sources extracted with 
a typical detection cell a few arc sec in size, the 
large extent of the galaxy and the faintness of the 
diffuse emission require a careful analysis of the 
background. The D25 circle* is ~14 arcmin in size 
(Tully 1988) and is larger than the ACIS S3 CCD. 
The large extent of the galaxy prevents the use of 
portions of the S3 chip outside of the galaxy as a 
measure of the background. We instead obtained 
the background from a "blank field" observation^ . 
The data were filtered identically by energy, then 
re-projected onto the sky using the aspect solution 
of the actual NGC 6946 observation. The expo- 
sure time of the background observation was ~300 
ksec. 

Point sources were detected as described in Holt 
et al. (2003), using a wavelet detection that is di- 
rectly tied to the point spread function's behavior 
with off-axis angle (Freeman et al. 2002). The 
point source detection reached a limiting luminos- 
ity Lx ~3xl0^'' erg s'^ in the 0.3-5 keV band 
which corresponded to ^10 counts. Holt et al. 
(2003) deemed the source detection complete to 
~lxl0^^ erg s~^; the completeness limit is lower 
by ^30% near the center of the CCD because of 
the sharper point spread function there where ~7- 
8 counts constitutes a detection. Most of the de- 
tected diffuse emission lies within 2.5-3 arcmin of 
the aimpoint. We do not see a trend of increasing 
numbers of weak sources as we approach the aim- 
point, suggesting a relative lack of contamination 
of the diffuse emission by weak point sources. We 

^The D26 circle is the diameter of the optical isophote at a 

surface brightness of 25 mag/arcsec^. 
®A description as well as links to the relevant files may be 

found at http://cxc.harvard.edu/contrib/maxim/ax;isbg/. 
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postpone until §5 a longer discnssion of the impact 
of possible unresolved point sources. 

Sources were removed from the data using the 
CIAO tool dmf ilth which cuts out each source, 
then replaces the resulting 'hole' by sampling the 
events in annuli surrounding the locations of the 
sources. The extraction apertures used enclosed 
>95% of the encircled energy. The annuli had in- 
ner radii 25% larger than the radii of the source 
extraction circles. Care was taken to ensure that 
the annular sampling did not include any source 
events. Figure 1 shows the results of all the steps 
just described: the energy-filtered data, minus 
point sources, and binned into 4-arcsec pixels. Im- 
mediately from the raw data (Figure 1) we see the 
diffuse emission as well as bright areas, particu- 
larly on the N and W arms and surrounding the 
nucleus. 

We then adaptively smoothed the data using 
the CIAO code csmooth. We adopted a gaussian 
smoothing kernel and selected the smoothing pa- 
rameter to be 3(7 as calculated locally. Finally, 

we subtracted a constant from the result to re- 
move the overall background as determined by the 
'blank field' background. Figure 2 shows the re- 
sults in the three bands. The 0.3-0.6 keV band 
shows little emission specific to NGC 6946 because 
of the high column in the galaxy's direction. The 
0.6-0.9 and 0.9-2 keV bands appear similar: con- 
tours fall in similar regions of the galaxy in each 
band. In the 0.9-2 keV image, the apparent emis- 
sion that extends beyond the region defined by the 
0.6-0.9 keV image is an artifact of the adaptive 
smoothing, driven by the lower counts in the 0.9-2 
keV band. We henceforth use just the 0.6-0.9 keV 
band for any comparisons. The "point sources" 
apparent in this figure do not correspond with any 
of the detected sources listed in Holt et al. (2003), 
nor do they correspond with point sources lying 
below the threshold. An examination of the raw 
data at full resolution shows a broad distribution 
of counts suggestive of a diffuse but clumpy source. 

3. Multiwavelength Comparison 

In this section we compare the X-ray image of 
the diffuse emission with images obtained in other 
bandpasses. 

The apparent "point sources" discussed at the 
end of the previous section do not correspond to 



any of the optically-identified supernova remnants 
(Matonick & Fesen 1997). The regions may corre- 
spond with prominent H II regions identified in 
the submillimeter band by Alton et al. (2002). 
The extracted counts from these regions yield lu- 
minosities of ~1.5xlO^'' erg using a simple 
thermal bremsstrahlung model with temperature 
of 0.3 keV (§4.4). 

Figure 3 shows the contours of the 0.6-0.9 keV 
emission overlaid on an Ha image. Note that 
in general the X-rays follow the Ha emission. 
This behavior is most easily visible on the promi- 
nent northern arm or at the sharp south-southeast 
edge. The weak arm to the northwest is not de- 
tected; nor is the weaker portion of the inner 
eastern arm. The matching behavior means that 
the X-ray emission arises from a hot diffuse ISM 
component or from the summed emission of point 
sources confined to the arms. We will return to 
this question in the discussion section. 

Figure 4 shows the 0.6-0.9 keV X-ray contours 
over the 6-cm polarized intensity map from (Beck 
& Hoernes 1996) and discussed in Frick et al. 
(2001). Beck & Hoernes (1996) showed that the 
regular magnetic field is strongest in the interarm 
region. In the previous figure, we established that 
the X-ray emission follows the northern arm. As 
in Beck & Hoernes (1996), it is clear the polarized 
emission lies external to that arm, so the X-ray 
image is consistent with their result. 

Figure 5 shows the 0.6-0.9 keV X-ray contours 
overlaid on the 850 fim map of Alton et al. (2002) . 
At 850 jdm, the nucleus as well as 3 knots of emis- 
sion are prominent; the knots lie to the northeast, 
east, and south-southeast of the nucleus. Overall, 
the regions of strongest X-ray emission correspond 
with the regions brightest at 850 /xm. These knots 
correspond to prominent H II regions. The 850 
/xm emission largely arises from dust heated by 
starlight and thermally radiating (e.g., Israel et 
al. 1999). 

Figure 6 displays the 0.6-0.9 keV X-ray con- 
tours on the 3 cm thermal radio emission (avail- 
able in Frick et al. 2001). In general, knots of 
X-ray emission align with the radio knots, partic- 
ularly on the NE arm. A longer X-ray exposure is 
required to ascertain whether the outer regions in 
the 3 cm image have a corresponding X-ray coun- 
terpart. 
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4. Spectral Fit: Introduction 

For all of the analysis, we fit the spectra of 
the diffuse emission + background and the back- 
ground simultaneously. This introduction de- 
scribes the extraction of the spectra and a fit to 
the background. 

Two approaches are generally used to extract 
a spectrum spread across a large expanse of the 
ACIS chips: (i) extract the data as small arrays, 
build a response matrix specific to each array, and 
sum the arrays; or (ii) correct the entire chip for 
the charge transfer inefficiencies (CTI), extract 
a single source-|-background spectrum, and con- 
struct one response matrix (the Penn State ap- 
proach, using the CTI corrector as described in 
Townsley et al. 2002). We adopted the second 
method. A correction for the time-dependent con- 
tamination was included by reducing the effective 
area at each affected energy^. 

Following the removal of the point sources, the 
spectrum of the diffuse cmission+background was 
extracted including all pixels within the S3 chip 
and with exposure >0.9 of the on-axis time (this 
eliminates edge effects). A total of 28248 counts 
were extracted representing the spectrum of the 
source plus background. A background spectrum 
was extracted from the reprojected blank field ob- 
servation using the same (x, y) image pixels as for 
NGC 6946 itself. The background spectrum con- 
tained 331600 coimts and 22334 counts after scal- 
ing the exposure time to that of NGC 6946. The 
source spectrum was adaptively binned so that the 
resulting spectrum contained at least 20 counts 
per channel. 

Figure 7 shows the two extracted spectra as 
an overlay. From this figure, it is clear that the 
galaxy's diffuse emission fades into the background 
near 2 keV and is blocked below ^0.5 keV by the 
column toward NGC 6946. This figure reinforces 
one's visual impression as well as quantitative ex- 
pectation: most of the visually apparent emission 
appears in the 1-2 kcV band where Chandra's ef- 
fective area takes on the largest values. The visi- 
ble 'emission lines' arise from particle background 
stimulation of material in the detectors; the large 
line at 1.8 keV is Si Ka (CXC Proposers' Guide 

'Information about the QE degradation is available at 
http://asc.haj:vaj:d.edu/cal/Acis/Cal_prods/qeDeg/index.html. 



2002)^. The difference between the two spectra of 
Figure 7 represents the net diffuse spectrum and 
contains 5914 counts. 

Figure 8 shows a fit to the background using 
a number of power law and gaussian components. 
The residuals are flat to within ~15% across the 
bandpass. Wc use this background model in all 
subsequent spectral fitting, freezing the parame- 
ters at their best-fit values. We justify this ap- 
proach by noting the background spectrum may 
be fitted to very high precision, it having been ac- 
cumulated from numerous long exposures. 

4.1. Fit to All Diffuse Emission 

Simple continuum models did not provide a 
good fit and in all cases left clear residuals. For 
example, the use of a single Mckal model pro- 
duced the lowest x^/^ of 1-21 of all the single- 
component models, but left clear residuals in the 
0.9 keV band. The residuals essentially force con- 
sideration of additional degrees of freedom. 

Four models provided essentially comparable 
goodness-of-fit values. At the interpretation stage, 
we will use the results from one or more models 
because any one of them could be deemed "cor- 
rect." First, we adopted a double 'Mekal' spectral 
model representing thermal emission from a dif- 
fuse gas as calculated by Mewc, Gronenschild, & 
van den Oord (1985); Mewc, Lemen, & van den 
Oord (1986) and Kaastra (1992) with Fe L calcu- 
lations supplied by Liedahl, Osterheld, & Gold- 
stein (1995). Second (and third), we used sin- 
gle and dual variable abundance Mekal models, 
allowing the abundances of species with lines in 
the 0.5-2 keV band to vary. Last, we fit the spec- 
trum with a simple thermal bremsstrahlung to de- 
fine the continuum and included as many gaus- 
sian lines as required by the data. For the multi- 
gaussian approach, the choice of continuum model 
is largely irrelevant, as any simple model will pro- 
duce a fit. A power law, for example, works as 
well as a bremsstrahlung model (the resulting in- 
dex is a steep ~7±0.55, 90% error); we adopt the 
bremsstrahlung model to provide a check on the 
fitted temperatures from the more complex mod- 

'^The online version shows the background spectrum in red, 
the sourcc-|-background spectrum in black; from the fig- 
ure, it is clear the source spectrum contains excess 1.8 keV 
emission above the instrumental feature 
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els. For the gaussian lines, we set their widths to 
zero, because the lines are unresolved with ACIS 
but broadened by the instrumental response, and 
their centers to the positions of expected lines 
(§4.3). The multi-gaussian approach provides in- 
formation on the ionization states contributing to 
the total emission. Figure 9 (top) shows the dual 
variable Mekal fit while the bottom figure shows 
the multi-gauss fit. 

We also included a power law component to ac- 
count for any excess emission between the source 
plus background and background spectra. The 
residual flux from this power law 'background' 
component is <6% of the total diffuse flux, veri- 
fying the quality of the background definition and 
its separation from the diffuse emission of NGC 
6946. Table 1 lists the parameters for the fitted 
models. 

4.2. Continuum Emission 

The continuum temperatures of the single con- 
tinuum models or the low temperature component 
of the dual models are identical within the errors 
with kT ~0.25-0.32 keV. Figure 10 (top) shows 
the variable Mekal model contours for the low- 
temperature component. The fitted column den- 
sity is also identical within the errors for all models 
with a range of 3.6-4.1 x 10^^ cm~^. These columns 
lie about a factor of 2 above the Dickey & Lockman 
column (~2xl0^^ cm~^. Dickey & Lockman 1990) 
but within the errors of the Schlegel, Finkbeiner, 
& Davis (1998) value of ~3-5xl02i cm-^. 

The integrated flux in the 0.5-2 keV band 
is ~3.7xl0~^^ ergs s~^ cm~^ (absorbed) or 
~2.5xl0~^^ (unabsorbed) . The integrated lu- 
minosity values are ^1.5x10^^ and ~lxl0^° erg 
s"^, respectively, or '~2.9xl0^'^ and ^1.9x10^* 
erg s~^ arcmin"^, respectively. 

These values differ somewhat from previous 
measurements. The ROSAT measure of the col- 
umn was low by about a factor of 2 (Schlegel 1994). 
The fitted temperature is also lower than the 
ROSATPSPC value of 0.55i^;25 keV; even though 
the error bars just overlap, the higher temperature 
behavior of the ROSAT fit is not unexpected. The 
fitted temperature, for example, from the PSPC 
data was almost certainly artificially high because 
it included photons from sources with harder spec- 
tra, such as X-ray binaries, that were incompletely 



removed. In addition, §4.1 shows that two spec- 
tral components are required for a good fit with 
one at kT ~0.7 while the PSPC data were consis- 
tent with a single temperature because of the de- 
tector's lower spectral resolution. The integrated 
flux is higher than the ROSAT value by ~25%, 
as is expected given the higher value for the col- 
umn density. As an aside, the very broad PSF of 
the ASCA mirrors prohibited an analysis of the 
diffuse emission from the NGC 6946 observation 
(Schlegel, Blair, & Fesen 2000). 

4.3. Line Emission and Abundances 

The abundances of the single variable Mekal 
model were first varied individually to test the size 
of the resulting error. An abundance, for which 
the 1 (7 error was consistent with 1.0, was reset to 
and fixed at 1.0. Two abundances (Si, Fe) were 
found to diff'er significantly from 1 (Figure 11). 
For Fe, the abundance is 0.67±0.13 solar at the 
99% level. For Si, the abundance is 3.0 with the 
90% contour just above 1.0 at 1.05; at 99%, the 
Si abundance is consistent with solar. No flares 
occurred during the observation, so the line is un- 
likely instrumental in origin. We return to this 
point momentarily. 

The gaussian centers were flxed at the ener- 
gies of potentially prominent lines, corresponding 
to O VII 0.57 keV, O VIII 0.653 keV, Fe I (L) 
0.705 keV, Fe XVII 0.826 keV, Ne IX 0.922 keV, 
Fe XX 0.996 keV, Ne X 1.022 keV, and Mg XI 1.34 
keV. Lines at energies lower than ~0.5 keV were 
deemed irrelevant because of the high column to- 
ward NGC 6946. Several lines were significantly 
detected at >99% (Table 1). At the energies of 
these lines, the ACIS spectral resolution is ~110- 
120 eV (CXC Proposers' Guide 2002). 

While we fixed the line centers at specific ener- 
gies to test for the presence or absence of known 
lines, we relaxed that requirement for the Si line to 
test whether the feature in the source-l-background 
spectrum was consistent with the background fea- 
ture. The background line^ is best fit with a gaus- 
sian at 1.775±0.002 keV (blend of Si Ka) while the 
diffuse Si line has a best-fit energy of 1.860±0.027 
keV (= Si XIII). 



Recall that the background spectrum may be arbitrarily 
precise because it is an accumulation of an arbitrary num- 
ber of blank sky events. 
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Assume for the moment the spectral fits are ac- 
curate measures of the hne emission. The presence 
of O VII and O VIII then provides a measure of the 
overall ionization in the diffuse emission. O VIII 
0.653 keV is the Lya transition; the blended O VII 
0.57 line(s) represent the ls2p to Is^ transition. 
The log of the ratio is -0.2±0.7 and has a large er- 
ror because of the decreasing continuum in the 0.5- 
0.6 keV band and because of the energy resolution 
of the CCD. Nevertheless, the ratio translates to 
limits of 9.6 < log Ugt <12.0, where Ug represents 
the electron density and t the cooling time (see, for 
example, Veddcr ct al. 1986). If the low 'vmekal' 
continuum temperature approximately represents 
the electron temperature, then it plus the line ra- 
tio restrict the O VII resonance/forbidden ratio 
to <0.6-0.8, a prediction for future observations 
using instruments with higher energy resolution. 
The values arc just barely consistent with ioniza- 
tion equilibrium and suggest non-equilibrium con- 
ditions exist in the diffuse medium of NGC 6946 
particularly given the lack of a good fit by the var- 
ious equilibrium ionization models. For the just- 
described electron temperature and line ratio val- 
ues, equilibrium requires ngt > 12.5-13.0. Lower 
values for the electron temperature plus lower val- 
ues for the interpreted line ratio push the esti- 
mated resonance/forbidden ratio toward higher 
values and increasing equilibrium conditions; for a 
given line ratio, higher electron temperatures in- 
dicate increasing non-equilibrium conditions. 

A recent paper by Page et al. (2003) presents 
the XMM-Newton RGS spectrum of gas in M81. 
We note that the spectrum shows the same re- 
solved oxygen emission lines with a ratio of -^0.65, 
indicating non-equilibrium conditions exist in the 
ISM of M81. While we acknowledge that the dif- 
fuse emission in NGC 6946 requires study with an 
instrument possessing higher spectral resolution at 
a spatial resolution comparable to Chandraa to 
confirm our analysis, the similarity between the 
line ratio in the RGS spectrum and our upper limit 
argues for non-equilibrium conditions in the ISM 
of both galaxies. 

A similar analysis using Ne IX and Ne X can 
not be carried out because the neon emission is 
blended with Fe XIX and Fe XX. The M81 spec- 
trum of Page et al. in fact illustrates that the 
blend occurs even at the resolution of the RGS. 

A deeper X-ray exposure will be necessary to 



provide sufficient events to tighten the line ratio 
and continuum temperature values to strengthen 
the non-equilibrium suggestion. A future mission 
with high spatial and spectral resolution is neces- 
sary to demonstrate conclusively whether the dif- 
fuse emission exists in a non-equilibrium condition 
by measuring the resonance/forbidden ratio of the 
O VII lines. The expected large point spread func- 
tion of Astro- E2 will blur too many point sources 
into the diffuse emission to provide such a good 
measure. 

4.4. Diffuse Emission without the Knots 
and the Spectrum of the Knots 

As an experiment, we extracted a spectrum of 
the diffuse emission in the knots visible in the 
smoothed images as well as a spectrum of the dif- 
fuse emission minus the knot emission. Figure 12 
shows the results. The spectrum of the 'no knot' 
diffuse emission is essentially identical to the com- 
plete spectrum discussed previously. This is not 
particularly surprising given the relatively bright 
emission distributed over the face of the galaxy 
and the relatively small number of events in the 
knots. 

The lower portion of the figure shows the knot 
emission itself. We fit the spectrum to obtain the 
column and continuum temperature with a single 
variable Mekal model, then fit for line; abundances 
as well as using gaussians to represent emission 
lines (as shown in the figure). Table 2 lists the 
model parameters. The background forms an in- 
consequential portion of the flux given the small 
spatial coverage of the knots. The abundances 
show enhanced Si but depleted O, Ne, and Fe. 
Gaussian line centers correspond to O, Si, and Fe; 
the fourth line at ~0.95 keV is most likely a blend 
of Ne IX and Fe XX. The integrated flux of the 
knot emission is ~3x 10"^'^ erg s~^ cm~^ or about 
5-10% of the original diffuse emission. 

5. Discussion 

The spectra of the point sources with sufficient 
counts to yield good spectral fits, x^/j^ ^1-3 for 
f ^50, arc uniformly harder than the diffuse spec- 
trum presented here (Holt et al. 2003). The lowest 
temperature of the best-fit bremsstrahlung mod- 
els for the point sources is ~1.9 keV. Similarly, 
the best-fit power law indices are typically ~2-3; 
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the steepest value is 4.9±0.4 versus the power law 
index fit to the diffuse emission of 7±0.55 (§4.1). 
These results all support the presence of a distinct 

spectral component distributed across the galaxy. 

Schlegel (1994) argued that the detection of 
the diffuse emission represented the hot ISM com- 
ponent and not the summed emission of point 
sources based upon a plausibility argument from 
the known luminosities of typical constituents. 
The argument is strengthened here because of the 
improved separation of point sources and diffuse 
emission as well as the deeper probe of the lu- 
minosity distribution of the Chandra observation. 
Holt et al. (2003) detected 72 point sources to 
~7xl0^^ erg s~^ and judged it complete to ~10^^ 
erg which is a factor of ^10 deeper than the 
PSPC data and a factor of ^7 deeper than the 
HRI image (Schlegel 1994; Schlegel, Blair, & Fe- 
sen 2000). 

Our argument can be strengthened still further 
by extrapolating the luminosity distribution to an 
arbitrary, low value of Lx. We adopt the lumi- 
nosity distribution of Holt et al. but eliminate 
from that distribution any sources not detected in 
the 0.5-2.0 keV band. A total of 7 sources are re- 
moved. A fit to the resulting distribution changes 
the distribution's slope marginally from 0.64±0.02 
to 0.62±0.01. If we extrapolate this distribution 
to Lx ~10^^ erg s""'^, we obtain a summed Lx 
of ~3.3xl0^* erg s""'^; pushing the distribution fit 
parameters to their positive error limits yields a 
summed Lx of ~3.6xl0^* erg s~^. When sub- 
tracted from the fitted diffuse luminosity, we ob- 
tain the lower bound on the diffuse luminosity, 
~l.lxl0^^ erg s'^. Of the events described in this 
paper as 'diffuse' thou, up to ^25% could originate 
in unresolved point sources below the detection 
threshold, leaving ^75% for the hot component of 
the ISM of NGC 6946. 

The above discussion assumes the lack of a 
yet-steeper contribution from the unresolved point 
sources. We can not assess the impact of such a 
component but note in defense of our argument 
that the luminosity distribution of the Milky Way 
does not contain evidence for the existence of such 
a component (Grimm et al. 2002). 

The unresolved point source contribution may 
very well be less than our 25% estimate. Grimm 
et al. (2002) estimate the number of sources above 
lO^'' erg in the Milky Way as -^700 with a pos- 



sible factor of 2 uncertainty. For NGC 6946, our 
extrapolation leads to ~3800 sources above 10^** 
erg s~^, a factor of 5 greater than the Milky Way. 
If the extrapolation for NGC 6946 lies closer to the 
Milky Way observations, the unresolved point con- 
tribution likewise drops, perhaps by as much as a 
factor of 2. Of course, this argument could equally 
well be pushed in the opposite direction by a sim- 
ilar factor because NGC 6946 is not a clone of the 
Milky Way. For example, NGC 6946 is extremely 
CO-rich, radio-bright, and undergoing very active 
star formation (Casioli et al. 1990; Walsh et al. 
2002). If we doubled the unresolved point source 
contribution, our estimate for the hot ISM contri- 
bution decreases to ~50% of the X-rays presented 
here as the 'diffuse spectrum' for a luminosity of 
~7xl03* erg s~^. 

Noting the caveats above but taking the original 
measurement at face value, our results represent a 
solid detection of the hot diffuse component from 
NGC 6946 free of any concerns over the resolved 
sources' point spread function wings that plagued 
results from ROSAT, for example (Schlegel 1994). 
Summing the unresolved point component, the dif- 
fuse component (with a luminosity in the range of 
'~50-75% of the diffuse Lx discussed in this paper), 
and the resolved sources (Holt et al. 2003) yields 
a total Lx of ~1.3xl0^° erg s'^ in the 0.5-2 keV 
band. The hot diffuse component then represents 
at least -^5-8% of the total X-ray emission of NGC 
6946 and potentially as much as ~10%. 

If we approximate NGC 6946, based on the 0.6- 
0.9 keV image, as a thin rectangle of length '^9'.5, 
width ^4'. 5, and thickness hkpcXl kpc, the result- 
ing volume V is ^3.8x10^^ hkpc cm'^. We adopt 1 
kpc for the unit of thickness for convenience. The 
X-ray-emitting gas may be distributed in a layer 
of different thickness. For example, Ehle & Beck 
(1993) found a consistent model for thermal radio 
emission and Faraday rotation with a thin disk, 
but argued for a thicker disk. They did so based 
upon the assumption that Faraday rotation and 
depolarization occurred in the identical volume, 
assumptions not necessarily valid^ because edge- 
on galaxies show both thin and thick disks in the 
thermal gas and synchrotron emission. 

The continuum model normalization parame- 

^Wc thank the referee for drawing our attention to this 
point. 
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ter is proportional to the emission measure £ = 
/rigdV. From the measured spectral tempera- 
ture (the low temperature component) and the 
emission measure, we may estimate additional pa- 
rameters describing the ISM following Wang et 
al. (1995) or Summers et al. (2003) to be consis- 
tent with their definitions. For a filling factor f 
and proton mass mp, the estimated gas density 
nl=£ I (fVhkpc) or ne-0.012 (fhkpc)"' cm'^. 
The mass of the gas in the diffuse component 
is M ~ fue V mp ~3.8xl0^ (hkpcf)5M0. The 
gas pressure p = 2 ne kT = 9.1x10"^^ (fhkpc)"^ 
dyn cm~^, the thermal energy E ^3 Ug kT V 
~ 5.2x10^^ (ii^^)5 erg, and the cooling time is 

tcooi ~3kT/Ane -1.9x10^ (fhkpc)^ yr, with A = 
1.6xl0~^^ erg cm^ adopted for the gas emis- 
sivity (Sutherland & Dopita 1993). The mass of 
X-ray emitting gas represents ~ 0.3^^^%, for a 
filling factor f, of the total molecular gas mass (H2 
~10i° Mq, Young & ScoviUe 1982). The cool- 
ing rate Mcooi ~ «^ (Nulsen, Stewart, & 
Fabian 1984) and for the measured values yields 
^2 Mr;;, yr^""^. If we assume the bremsstrahhmg 
or Mekal temperatures represent the actual tem- 
peratures of hot gas, they correspond to ther- 
mal velocities of 180-300 km s^^. The number 
density is about the same, the thermal energy is 
-^5 times lower, the pressure ~20% lower, the X- 
ray-emitting gas mass ^6 times higher, and the 
cooling time about 2 times shorter, than the cor- 
responding values for NGC 4449 based upon its 
Chandra observation (Summers et al. 2003). How 
these values correlate with other galaxian parame- 
ters requires a considerably broader survey of dif- 
fuse emission. At this point, we note that NGC 
6946 has a higher total mass, more H I gas, and a 
higher 60 /itm luminosity (Young et al. 1989), in- 
dicating more vigorous star formation, than NGC 
4449 possesses. 

The analysis of the RGS spectrum of M81 (Page 
et al. 2003) finds a dual temperature plasma^° 
describes most of the spectrum with kTi = 
0.18±0.04 and kT2=0.64±0.04 keV, very simi- 
lar to our results achieved at a lower spectral but 
higher spatial resolution. They infer a cooling 
time of ~4xl0^ ^/f yr, similar to our inferred 

'Page et al. formally require a contribution from a third 
component with kTs = l.T^g'g keV. They argue that this 
component arises for the most part in the bulge of M81. 



value. Kuntz et al. (2003) obtain a best-fit dual 
temperature model of kT=0.20 and 0.75 keV for 
the diffuse component in MlOl. Finally, Ehle et 
al. (1998) fit the diffuse component of M83 with 
a dual temperature model and obtained kTi ~0.2 
and kT2 ~0.5 keV. These results all support the 
existence of at least a 2-temperature hot ISM in 
spiral galaxies. 

The X-ray knots correspond spatially with local 
maxima in the 850 fim image. The 850 /im max- 
ima are knots of dust heated by starlight and cor- 
respond with H II regions or areas of star forma- 
tion activity (Alton et al. 2002). That the X-ray 
image shows local maxima in the same location 
raises the possibility that the knots are blowout 
regions (e.g. Ott et al. 2001). Walsh et al. (2002) 
showed that the total 6 cm radio image correlated 
with images from the mid-IR and argued that the 
apparent single correlation was in fact two corre- 
lations, between the clumpy mid-IR vs. 6 cm ther- 
mal radio and the diffuse mid-IR vs. 6 cm non- 
thermal radio emission. The diffuse radio emission 
is dominated by the nonthermal component; the 
radio emission from clumps was largely thermal in 
origin. The correlation with the X-ray then sug- 
gests the X-ray emission is thermal in origin, so 
line emission should be present. The knot emis- 
sion shows at best weak evidence for line emis- 
sion within the statistical uncertainties. A higher 
signal-to-noise spectrum is certainly needed to im- 
prove or refute the upper limits on the presence of 
X-ray line emission. We intend to pursue addi- 
tional observations, particularly deeper observa- 
tions. 

The authors thank the anonymous referee for 
suggestions that improved the presentation of the 
paper. This research was supported by contract 
number NAS8-39073 to SAO. The work made use 
of the CI AO software built under contract to SAO. 
The authors thank W. Walsh for providing the CO 
data, R. Beck for the 3, 6, and 6 cm PI data, R. 
Fesen for the Ha image, and F. Combes for the 
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Table 1: Fitted Parameters for the Spectrum of the Diffuse Emission** 



Modcl^' 



Parameter N 



H 



U.ZOO_o.o25 

' '-'-0.08 
2 -10+0.48 



+0.10 
-0.08 



Norm Fhix'tr EqW X^/'^ 



Mckal (kXi) + 
Mekal (kXa) + 
Power Law (Bkg) (F) 



Brems (kT) + 
gaussian (norm, 
gaussian (norm, 
gaussian (norm, 
gaussian (norm, 
gaussian (norm, 
gaussian (norm, 
gaussian (norm, 
gaussian (norm, 
gaussian (norm, 
Power Law (Bkg 



O VII) + 

O VIII) + 
Fe I L) + 
Fe XVII) + 

Nc IX) + 
Fe XX) + 
NeX) + 
Mg XI) + 
Si XIII) + 

) (r) 



VarMekal (kTi) + 
VarMekal (kTa) + 
abundance: Si + 
Power Law (Bkg) (F) 

SingVarMekal (kT) + 
abundance: O + 
abundance: Si + 
abundance: Fc + 
Power Law (Bkg) (F) 



+0.30 
-0.12 



0.32 
0.569 

0.653 
0.705 
0.826 
0.922 
0.996 
1.022 
1.340 
1.860 
2.24^ ' 

0-245_o.„i5 
0_7^+o.io 



0.25 



-0.050 



-0.08 
-1.44 
-1.19 
9 09+0-53 
^•'-'^-0.32 



2.26 



0.276tH^^ 
0.83t°;?° 

o 00+1-95 
^■"^-1.90 

67+° " 

'-'•O' -0.12 
9 09+0-33 
-^•'-'^-0.25 



0.36 



0.41 



+0.10 
-0.12 



0.40 



+0.07 
-0.17 



9.4C-4 
1.3e-4 
4.3e-5 

2.3e-2 



4.5e-5 

1.3e-3 
l.le-4 



2.3±0.4c-12 
1.8±0.6e-13 
1.0±0.3e-13 

2.3±0.5e-12 
403+90" 

'±UO_220 
^""-203 

<600 

or 0+90 

000+168 
^°^-120 

<100 

QO+80 



2501^^° 
9.5±0.4e-14 



2.5±0.6e-12 
2.8±0.8e-13 



4.2e-5 9.4±1.9e-14 



0.39^^5 l-5e-3 2.4±0.6e-12 1 



5.5e-5 1.2±0.2e-13 



22 351 



14 342 



19 351 



23 350 



"Errors are all 90%. Upper limits are all 99%. The 'Bkg' power law is used to account for any residual background emission 
incorrectly subtracted. 

'All gaussian line widths were fixed to 0.0 because the lines are unresolved but broadened by the instrumental resolution; line 
centers are fixed at the listed values. 
=Units = 10^2 cm-2. 

''Unabsorbed flux in 0.5-2.0 keV band with exponent indicated by e-XX = 10"'''^''^; units = erg cm"^; EqW = Equivalent 
Width, units = eV. 

''The degrees of freedom represent the total number in the source spectrum; the background degrees of freedom, 558, have been 
subtracted from the total number of degrees of freedom counted by the software. 
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Table 2: Fitted Parameters for the Knots/No Knots Spectra** 
Model'' Parameter Nh'' 



Norm Fliix'tr EqW -y^jv 



Knots Emission: 
Brems + 

gaussian (O) + 0.652 

gaussian (Fe) + 0.822 
gaussian (Ne+Fe) + 0.95 
gaussian (Si) + 1.80 

VMekal + 
abundance: O + 
abundance: Nc + 
abundance: Si + 
abundance: Fe + 



20 



+0.08 
-0.02 



24+° °^ 
0.37t°-i^ 
0.39t°i^ 
2.86t?;?^ 



fi2+° °^ 

U.OZ_Q J^Q 

l.le-4 

2.7e-5 
7.4e-6 
7.5e-7 



0.40 



+0.07 
-0.17 



4.1e-3 
180±110 

170±75 

90±45 

<500 

2.2e-5 



0±1.8e-13 



0.75 310 



5±0.9e-13 0.76 315 



No Knots Emission: 
Brems + 

gaussian (norm, O VIII ) + 
gaussian (norm, Fe I L ) + 
gaussian (norm, Fe XVII ) + 
gaussian (norm, Si XIII ) + 
Power Law (Bkg) 



0.28t°;°^ 

0.653 

0.707 

0.826 

1.860 



2.09 



+0.36 
-0.27 



tiO+o O"* 

U.OU_Q gg 



7 


05e-3 


3.6±0.94e-12 


3 


95e-4 


265±145 


3 


69e-4 


245±85 


2 


02e-4 


215±70 


5 


33e-6 


440±170 


3 


4e-4 


2.4±l.le-13 



25 347 



VMekal + 
abundance: Ne + 
abundance: Si + 
abundance: Fc + 
Power Law (Bkg) 



0.28^°:^^ 
0.74t°;?i 

o qc + 1.95 

fS3+°-^^ 

^■O-J-O-QO 
^•'-'^-0.25 



n 40+0.05 



2.57e-3 4.3±0.94e-12 1.30 351 



3.2e-4 2.4±l.le-13 



"Errors are all 90%. Upper limits are all 99%. The 'Bkg' power law is used to account for any residual background emission 
incorrectly subtracted. 

'AH gaussian line widths were fixed to 0.0 because the lines are unresolved but broadened by the instrumental resolution; line 

centers are fixed as listed. 
'^Units = 10^2 em"2^ 

''Unabsorbed flux in 0.5-2.0 keV band with exponent indicated by e-XX = 10"-''-'''-; units = erg s~^ cm~^; EqW = Equivalent 
Width, units = eV. 

'^The degrees of freedom represent the total number in the source spectrum; the background degrees of freedom, 558, have been 
subtracted from the total number of degrees of freedom counted by the software. 
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Fig. 1. — Chandra observation of NGC 6946 with detected sources removed. The data have binned by a 
factor of 8 to create 4 arcsec pixels. FIGURE IN SEPARATE JPG FILE: fl.jpg 



Table 3: Model Parameters for Instrumental Background 



Component 


Parameter*^ 


Norm 


Gauss (7 (eV) 


Power law 


-0.223 


3.70e- 


■5 




Power law 


3.871 


3.47c- 


■5 




Power law 


2.789 


6.77c- 


■6 




Gaussian 


0.670 


1.03c- 


■5 


4.84e-4 


Gaussian 


10.51 


9.12c- 


3 


1.642 


Gaussian 


2.149 


2.21c- 


■5 


5.64C-2 


Gaussian 


0.577 


4.09c- 


■5 


3.65e-2 


Gaussian 


0.232 


5.08c- 


3 


4.77e-2 


Gaussian 


1.775 


1.83c- 


■5 


2.61C-2 


Gaussian 


2.629 


2.75c- 


■5 


0.468 


Gaussian 


0.882 


4.06e- 


6 


5.06e-2 



"for power law, parameter = F; for gaussian, parameter = line energy 

Fig. 2. — Mosaic of the subbands of the diffuse emission, all adaptively smoothed with a gaussian kernel to 
place the minimum background S/N at 3.0. The 0.3-0.6 0.6-0.9, and 0.9-2.0 keV bands are shown contours 
for each band overlaid. The 0.3-0.6 keV band shows little emission because of the high column density; 
contour values are 0.34, 0.38, 45, and 0.55 counts s~^ pixcl"-'^. Contours in the 0.6-0.9 band lie at 0.45, 0.60, 
0.75, 0.90, 1.05, and 1.20 counts s'^ pixel'^ Contours in the 0.9-2.0 band lie at 0.50, 0.65, 0.80, 0.95, 1.20, 
and 1.60 counts s~^ pixel"^. The background levels are 0.26, 0.22, and 0.09 counts s~^ pixel"^, for the 3 
bands, respectively. The short vertical bar in the lower right represents a 2 arcmin scale bar. FIGURE 
IN SEPARATE JPG FILE: f2.jpg 



Fig. 3. — X-ray contours for 0.6-0.9 keV overlaid on an Ha image; the Hq image is courtesy of Robert Fesen 
(Matonick & Fesen 1997). The X-ray contours are described in Figure 2. FIGURE IN SEPARATE JPG 
FILE: f3.jpg 



Fig. 4. X-ray contours for 0.6-0.9 keV overlaid on 6 cm Polarized Intensity map from Beck & Hocrnes 
(1996). Note the contours lie between the 'polarized spiral arms' as described by Beck & Hoernes (1996). 
The vertical bar is 2' in length. FIGURE IN SEPARATE JPG FILE: f4.jpg 



Fig. 5. SCUBA 850 /i image with contours of X-ray diffuse emission (0.6-0.9 keV) overlaid. The 850 /im 
data arc from Alton et al. (2002). The gray levels have been set to match approximately the contours defined 
by Alton et al. in which the lowest visible gray level corresponds to their 3a isophote. The vertical bar is 2' 
in length. FIGURE IN SEPARATE JPG FILE: f5.jpg 
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Fig. 6. — 3 cm thermal radio emission with contours of X-ray diffuse emission (0.6-0.9 keV) overlaid. The 3 
cm data are from Frick et al. (2001). FIGURE IN SEPARATE JPG FILE: f6.jpg 



F 




0.5 1 2 

channel energy (keV) 



Fig. 7. — Plot of the diffuse -I- background (upper curve) and background (lower) spectra. The background 
spectrum was extracted from a blank sky observation after reprojecting those events using the aspect solution 
from the NGC 6946 observation and scaling by the ratio of the exposure times between the background and 
the NGC 6946 observations. The diffuse -I- background spectrum contains 28248 counts; the background 
spectrum contains 22334 counts. The difference is the net spectrum of the diffuse emission and contains 
5914 counts. 



Fig. 8. — Fitted background spectrum using a combination of power law and gaussian model components. 
This figure may be compared to the corresponding figure in CXC Proposers' Guide (2002). The residuals 
are flat to within -15%. FIGURE IN SEPARATE JPG FILE: fS.jpg 



Fig. 9. — (top) Spectral fit to the diffuse + background spectrum using a variable dual Mekal plasma model. 
The instrumental background has been fitted using the model of Figure 8. The individual model components 
are also shown, (bottom) Spectral fit to the diffuse data using a thermal bremsstrahlung and a series of 
gaussians. The fitted bremsstrahlung temperature is identical, within the errors, to the temperature of the 
top figure. FIGURES IN SEPARATE JPG FILES: f9a.jpg, f9b.jpg 
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NGC 6946; DualVMekal 



0.6 0.7 0.8 0.9 

high kT (keV) 

Fig. 10. — The contour plot for the variable Mekal models of the errors on the fitted parameters of the two 
thermal temperatures and the column. The contours for this and the next figure are 1 a, 90%, and 99% for 
2 parameters of interest (Ax^ = 2.30, 4.61, and 9.21). 



15 





Fig. 11. — Contours for the abundances of Fe and Si from the variable Mckal spectral fit. The abundance of 
Si differs from solar at the 90% level but not at 99%. The Fe abundance differs at 99%. 




0.5 1 2 

chonnel energy (keV) 

Fig. 12. — (upper) Spectral fit, using a thermal bremsstrahlung and a series of gaussians, to the diffuse data 

minus the knots of emission. The fitted bremsstrahlung temperature is identical to the bremsstrahlung tem- 
perature of figure 9 within the errors, (lower) Spectral fit to the emission of the knots using a bremsstrahlung 
and three gaussians. Because of the small spatial extent of the knots, the background is nearly negligible and 
has not been shown in the figure; it was included in the spectral fitting. FIGURE fl2a IN SEPARATE 
JPG FILE: fl2a.jpg 
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